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Employing a variational approach that takes into account electron-phonon 

and magnetic interactions in La\- x A x MnO^ perovskites with < x < 0.5, the 

effects of the magnetic field and the oxygen isotope substitution on the phase 

diagram, the electron-phonon correlation function and the infrared absorption 

at x = 0.3 are studied. The lattice displacements show a strong correlation 

with the conductivity and the magnetic properties of the system. Then the 

conductivity spectra are characterized by a marked sensitivity to the external 

parameters near the phase boundary. 
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The perovskite oxides La\- x A x MnO^ (A represents a divalent alkali element such as 
Sr or Ca) exhibit the colossal magnetoresistance (CMR) effect within the hole doping 
range 0.2 < x < 0.5. These compounds are characterized by a transition from a metal- 
lic ferromagnetic (FM) low-temperature phase to an insulating paramagnetic (PI) high- 
temperature phase that is associated with dramatic changes in their electronic and magnetic 
properties. 

These materials were first studied in the 1950's for their peculiar strong correlation be- 
tween magnetization and resistivity.! The ferromagnetic phase was explained by introducing 
the double exchange mechanism!'! that takes into account the combined effect of the Mn e g 
electron hopping between nearest neighbor sites and the very strong Hund's exchange with 
the localized Mn t 2s electron spins. Later on, in addition to the double-exchange physics, a 
strong interaction between electrons and lattice distortions has been suggested in order to 
explain the CMR phenomenon.! 

At high temperatures, above the Curie temperature T c , the Jahn- Teller polaron formation 
has been confirmed by many experimental measurements^ in particular by the activated be- 
havior of the conductivity,^!^ the thermopower,0'El the Hall mobility^ and more directly by 
atomic pair distribution,^ electron paramagnetic resonancejil x-ray and neutron scattering 
studiesiiUll X-ray-absorption fine-structure (XAFS) spectra have found significant changes 
in the local structure of these compounds showing the crossover from large to small polarons 
across the metal-insulator (MI) transition and the direct relationship between lattice distor- 
tions, charge distribution and magnetismJlH^l Indeed the charge carriers partly retain their 
polaronic character even below T c , as proved by neutron pair-distribution-function analysis0 
and resistivity measurements.^ Pseudogap featuresllzl and conductivity spectr ME have 
also been discussed in terms of a strong coupling to lattice distortions. Finally the existence 
of the strong electron-phonon coupling and polarons has been demonstrated by the giant 
isotope shift of the Curie temperature. Nowadays the key-role of the polaron formation 
in manganites is generally recognized. H@ 

Recently, studies have stressed the intrinsic strong tendencies toward the phase sepa- 



ration in manganites.iSHHi] Direct evidences for coexisting localized insulating and delo- 
calized metallic components have been reported from tunnelling spectroscopy, neutron and 
electron measurements.0'&0 The coexistence between hole-poor and hole-rich phases has 
been theoretically studied by using exact numerical approaches on small lattices assuming 
classical Jahn- Teller phonons.EH Within this approximation the optical propertiesSB have 
been studied but, at low temperatures, in the FM phase, the narrow Drude peak, that is 
experimentally seen,il'i! cannot be obtained. 

In a recent paper based on a variational approach,!^! some of us have shown that the 
quantum nature of the Jahn- Teller phonons and the polaron formation can be important to 
explain the experimentally observed tendency of manganites to form inhomogeneous mag- 
netic structures near the phase boundaries. Employing the scheme proposed in that work, in 
a subsequent paper@ spectral and optical properties of manganites have been deduced at dif- 
ferent temperatures for x = 0.3 (around x = 0.3 the CMR effect is very pronounced in many 
manganitesi ) . By using the formalism of generalized Matsubara Green's functionsE^^ the 
scattering rate of charges interacting with lattice and spin fluctuations has been determined 
and used to calculate the optical properties of the system. At low temperatures, in the FM 
phase the system is characterized by two types of optical response: a Drude term and a 
broad absorption band due respectively to the coherent and incoherent motion of large po- 
larons. At high T the infrared absorption is due to the incoherent small polaron dynamics. 
Upon cooling, the optical spectra have showed a transfer of spectral weight from high to low 
energies filling up the low frequency optical gap present in the PI phase in agreement with 
experimental data.!!©!! 

In this paper, within the proposed variational approach,@@ we calculate the lattice dis- 
placement probability distribution function (LDPDF), a quantity measurable by XAFS 
and introduced in order to study the lattice deformations due to the polaron formation.il 
It allows to estimate the Debye- Waller (DW) factor for several T at x = 0.3. Near and 
within the PI phase, the LDPDF shows that sites occupied by an electron (Mn 3+ ) are 
characterized by strong lattice distortions, while sites without conduction electrons (Mn 4+ ) 



by a vanishing average static displacement. Then we examine the effects of an external 
magnetic field on the phase diagram, the lattice distortions, the infrared conductivity and 
the resistivity of the system at x — 0.3. Near the phase boundary, the introduction of 
the field induces a transfer of spectral weight from high to low energies in the conductivity 
spectra and, consequently, large variations in the resistivity. In this approach, the CMR 
effect is due to the subtle balance between two coexisting phases characterized by different 
lattice, spin and electronic properties. The comparison between the DW factor in presence 
of the magnetic field and the magnetoresistance ratio points out the correlation between 
lattice distortions, conductivity and magnetic properties. Finally the isotopic substitution 
of the oxygen from 16 to 18 is faced, finding that it increases the electron-phonon cou- 
pling reducing the Curie temperature T c . In the vicinity of the MI transition, the isotope 
effect causes large changes of the conductivity, so that it represents another example of the 
sensitivity of the system to the variation in external parameters. 

In Sec. I and II the variational approach and the lattice distortions are respectively 
discussed; in Sec. Ill the calculation of the damping of the particle motion and the optical 
properties are reviewed,^ in Sec. IV and V the introduction of an external magnetic field 
and the isotope effect are respectively examined. 



I. VARIATIONAL APPROACH 

We adopt a model that takes into account the double-exchange mechanism, the coupling 
of the e g electrons to lattice distortions, the super-exchange interaction between neighboring 
localized t 2g electrons and the effects of an external magnetic field /z ea , t .@'EI The interaction 
to longitudinal optical phonons derives from the Jahn- Teller effect that splits the e g double 
degeneracy. Adopting the single orbital approximation (reasonable in the doping regime 
where CMR occurs@), the Hamiltonian reads 

H = -t ( %y^~f~ J ° ] i c i +UJ oJ2 a \ a i + 9UQ J2 c ld (<H + 4) 
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+e J2 & ' &3 ~ 9s^B J2 hext ■ Si- fJ,J2 c \°i- (!) 

<i,j> i i 

Here t is the transfer integral between nearest neighbor (nn) sites < i,j > for electrons 
occupying e g orbitals, S is the spin of the t 2g core states (S = 3/2), Sq 3 is the total spin 
of the subsystem consisting of the spins localized on a pair of nn sites and the conduction 
electron, cj (q) creates (destroys) an electron with spin parallel to the ionic spin at the i-th 
site. The first term of the Hamiltonian describes the double-exchange mechanism in the 
limit where the intra-atomic exchange integral J is far larger than the transfer integral t. 
In the second term of eq.(|I]) a\ (a») is the creation (annihilation) phonon operator at the 



site i, ujq = yk/M denotes the frequency of the optical phonon mode, with k restoring 
force per length unit of the local oscillator and M mass of an oxygen atom surrounding 
the manganese atom. The dimensionless parameter g indicates the strength of the electron- 
phonon interaction within the Holstein modelS 

- - 

where A is the energy per displacement unit due to the coupling of the charge carriers with 
the lattice. Finally in the Hamiltonian (HD e represents the antiferromagnetic super-exchange 
coupling between two nn t2 g spins, g s the Lande's electron spin factor (g s ~ 2), jj,B the Bohr 
magneton and u is the chemical potential. 

The hopping of electrons is supposed to take place between the equivalent nn sites of 
a simple cubic lattice separated by the distance |n — n'\ = a. The units are such that the 
Planck constant h — 1, the Boltzmann constant &b=1 and the lattice parameter a—1. 

We perform two successive canonical transformations to treat the electron-phonon inter- 
action variationally. The first is the variational Lang-Firsov unitary transformationil'il 



U\ = exp 



9 J2 (f c h + A ) ( a i - a ) 
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(3) 



where the parameter / measures the degree of the polaronic effects and the parameter 
A denotes a displacement field describing average static lattice distortions. The second 
Bogoliubov-type transformationil 



U 2 = exp 



j 



(4) 



with a variational parameter, introduces correlations between the emission of successive 
virtual phonons by the electrons. 

The transformed Hamiltonian H = U^UxHU^U^ * s difficult to treat exactly, so, in 
order to obtain a variational free energy,00 we introduce a test Hamiltonian characterized 
by electron, phonon and spin degrees of freedom non mutually interacting 

Htest = -t e ff c l c j +LjJ oYl <4°i + NuJ o sinh 2 (2a) + Nu g 2 A 2 

<i,j> i 

~9s^b J^h- S { + (r] - c Ui- (5) 

i i 

The quantity t e ff denotes the effective transfer integral 

where the symbol <> indicates a thermal average and the quantity St is 

St = g 2 f 2 e~' a (2iV + 1) , (7) 

with Nq the average number of phonons with frequency lDq = uj cosh(4a). In the test 
Hamiltonian (|5]), is the number of lattice sites, h the sum of the external magnetic field 
and the molecular magnetic field effective in a cell containing two neighboring sites^ and 
the quantity rj 

V = "o9 2 f (/ - 2) + 2u; g 2 (f - 1) A (8) 

measures the electronic band shift due to the electron-phonon interaction. 

We employ the Bogoliubov inequality in order to derive the variational free energy of the 
system 

F < F test + {H- H test ) t (9) 

where <>t indicates a thermodynamic average made using the test Hamiltonian, so that 
the local spin dynamics is studied within a variational mean field theoryill The free energy 
per site becomes 



^ = ftt + T log (l - e-**) + u sinh 2 (2a) + 1^(1 - /)V - T log ug 



±^ZS 2 m 2 s + TXM (10) 



where the electron free energy ff est is calculated considering the renormalized band = 
£k + Vj where the band dispersion is 

£k = — 2t e ff[cos(k x ) + cos(ky) + cos(k z )). (11) 

In eq.(|TDp (3 is the inverse of the temperature, v$ is the partition function of the localized 
spins, the top and bottom sign of e hold, respectively, for the ferromagnetic and antiferro- 
magnetic solutions of the localized spins and Z indicates the number of nearest neighbors. 
The quantity A is a dimensionless variational parameter proportional to the total magnetic 
field h 

A = /3g af i B Sh (12) 

and M = B(X) is the Brillouin function representing the ratio between the magnetization 
of the localized spins and the saturation magnetization. 

In order to simplify the calculations, we consider a semicircular density of states 

9^) = f-L) 0(W - \e\)VW^^ (13) 



where W = Zt e ff is the renormalized band half-width and 9{x) is the Heaviside function.!^ 2 "'' ! 



In the intermediate electron-phonon coupling regime, the free energy fllOf) gives rise to a 
region of coexisting phases characterized by different electron densities p\ and P2 and consti- 
tuted by large and small polarons, LP and SP, respectively. Thus, near the metal-insulator 
transition, the system segregates in ferromagnetic and antiferromagnetic or paramagnetic 
domains of LP and SP. 

The variational approach combined with the Maxwell construction allows to determine 
the fractions of volume of LP and SP. Within the regime of the macroscopic phase sep- 



aration (PS) ,E3E2l it is possible to calculate any thermodynamic property B of the system 



by means of the linear combination of the properties B\ and B2 characteristic of the single 
homogeneous phases assuming the respective volume fractions (Vi/V) and (V 2 /V) as weights 

B=(£W+f£K (14) 



V J \v 

Therefore we assume that in this regime the properties of the system are independent of 
morphology of coexisting domains neglecting the surface energy cost. In the field of CMR 
manganites, this two-fluid or two-component scenario has been proposed by other workers 
from a phenomenological point of viewlS Furthermore, adopting this assumption, in a 
previous work@ the results for spectral properties and infrared absorption have turned out 
in good agreement with experiments. 

In this paper the quantities are calculated taking into account the phase diagram^ 
evaluated by using the following values of the model parameters: t = 2u , g = 2.8 and 
e = 0.01*. 



II. LATTICE DISTORTIONS 

In this section we deal with the lattice distortions in absence of the external magnetic 
field. 

Within the above-mentioned variational approach, some of us have discussed the Jahn- 
Teller distortion of the Mn 3+ ionS It consists in an axial elongation of two MnO bonds of 
the MuOq octahedral Above a crossover temperature T* < T c , there are two distortions 
of this same type with a different degree of elongation along the axial direction, showing 
evidence of coexistence of large and small polarons. The average displacement of the site 
% from the equilibrium position, when one electron is on site i, is in agreement with the 
experimental data on the separations between local Mn — O bond lengths.^ 

Now we calculate the LDPDF P(X), that is more appropriate to make a comparison 
of the theory with the data obtained by XAFS measurements.il Considering the relevance 
of the axial elongation, we can confine ourselves to one spatial direction, so 
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1 f°° 

P(X) =< 5(X - X t ) >= — / dqe iqX < e- tqXl >, (15) 

where X^ = (a, + aj) / \/2Mwq is the displacement operator at the site i along a reference 
axis. It has been evaluated performing the two canonical transformations and making 
the decoupling^ in the electron and phonon terms through the introduction of H tes t ©• 
We get 

P {X) = -U \pe-^ 2 ' 12 + (1 - p)e-^ 2 ' l2 \ (16) 



where p is the electron density, I = /^?§^ b = 9\[w^ ( A + /) and c = # \/2? A - 
We obtain two terms: the first is due to the average occupation of the site i and involves 
directly the polaronic distortion controlled by the variational parameter /, the second is 
due to the remaining configurations without electron on the site % and it is related to the 
average lattice displacement A = p(l — /). In Fig. 1(a) we report the distribution P(X) 
for different temperatures at x — 0.3. In the regime of coexisting phases LDPDF takes 
contribution from both LP — FM and SP — PI distributions by means of the eq . (|T4]) . The 
LP — FM phase is characterized by a single-peak LDPDF centred about the average lattice 
displacement A (here / is smaller than A). With increasing temperature, the spreading of 
P{X) signals the PS regime and near T c two well-defined peaks appear giving a bimodal 
LDPDF . Near and in the insulating state, occupied sites (Mn 3+ ) have strong distortions 
while non-occupied sites (Mn 4+ ) are characterized by an average static displacement A 
equal to zero 

We have derived the mean-square deviation of this distribution Ax 2 =< x 2 > — < x > 2 
that provides the DW factor (see in Fig. 1(b) DW as a function of T at x = 0.3). The DW 
factor increases rapidly with temperature near T c in agreement with experimental data.0@0 



III. DAMPING AND OPTICAL PROPERTIES 



In this section we briefly review the calculations of the scattering rate and the infrared 
absorption due to lattice distortions and spin fluctuations.^ 
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Retaining only the dominant autocorrelation terms at the second step of iterationEH^], 
the self-energy £( 2 ) (k, ioj n ) and the scattering rate 

T(k) = f(k, u = 60 = -29fE?2 (k, u = 60 (17) 

are derived. Making the expansion into the series of multiphonon processes, the scattering 
rate is expressed as 

r(k) = r(60 = J- l—phon (60 + J- multi—phon 

(60 + r 5p j n _ FiMCt (^ k ) (is) 

where Y\- p h on is the contribution due to single phonon processes, T mu iu- P hcm represents the 
scattering rate by multiphonon processes and Ts P i n -Fiuct denotes the damping term by spin 
fluctuations. At low T, in the LP phase, single phonon emission and absorption represent 
the main mechanism of scattering. Furthermore, at higher temperatures the damping due 
to spin fluctuations is effective in the energy range around the chemical potential \x. 

The scattering rate turns out fundamental to derive the infrared absorption of the system. 
The real part of the conductivity is obtained by the current- current correlation function 

»*a» = ( i9) 

where the current operator j a suitable for the Hamiltonian ([[]) is 

j a = ite 8& 2 S + 1^ ') c i+ s & Ci - ( 20 ) 
The conductivity can be expressed as sum of two terms 
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n™^ 1 ) (lu) + n ret w (o, 

= - " — >l = ^ { « a a d) (") + ^T h \uj). (21) 



The first term 3?cr^ rid - ) represents the band conductivity because the absorption is not ac- 
companied by processes changing the number of phonons. On the other hand, the second 
term marked by the apex "incoherent" in eq. (|2T| ) derives from inelastic scattering processes 
of emission and absorption of phonons. The band conductivity is derived as 
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where C(£,u) is 



and h(£) reads 



1 



r(0 



r 2 (0+o; 2 



(23) 



(24) 



The latter term of the conductivity becomes 



'2e 2 t 2 ' 
'2e 2 t 2 ' 



-2S 7 



-2S T 



' S +1/2 " 
25+1. 



/+oo 
-co 



'5o + l/2 s 



'^o + l/2 s 



2S+1 J I \\ 2S + 1 ) I 

/+oo r+oo 
d£ / d^g(Og^i)[n F ^-u)-n F (0]C^,^,u) 
-CO J —CO 

where is the density of states ( |T3"D and the function i?(£,£i,u;) is given by 



fl(e,ei,w) = 2S/,(z)sinh {^f 



We notice that J/(£, £i, k>) 



(25) 



(26) 



J,(f , 6, w) = C(£, 6, w + too) M£ - too - u) - n F (£ - too)] [N (lco ) + n F (£)] (27) 
and u;) 

6, w) = C(£, Ci, - to )[n F (e + too - w) - n F (e + to )] [iV (to ) + 1 - n F (£)] (28) 

describe phonon absorption and emission processes, respectively. 

In the limit of high temperatures (T > 0.39co>o) and for SP excitations, the incoherent 
absorption is prevalent. In this case the conductivity consists of a sum of narrow Lorentzian 
functions centred on the points uluq respectively.il A derivation suitable for high tempera- 
tures can be performed and gives 



[ u J \[ 2S+1 



?{ 



exp 



Az 



exp 



ft I -\2 

Tz {uJ + z) 



p(l - p) (29) 
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where p is the electron concentration. 

The conductivity cr QjCf is given by eq.(pl|) and eq.([29|) in LP — FM and SP—PI phase, re- 
spectively. In the regime of coexisting phases, the two preceding conductivities are combined 
by means of the eq.([U|). 

At low temperatures, in the FM phase the system shows two types of optical response: 
a Drude term and a broad absorption band due respectively to the coherent and incoherent 
motion of large polarons. With increasing T, the optical spectra are characterized by a 
transfer of spectral weight from low to high energies. At high T the infrared absorption 
consists in a band peaked approximatively around the energy 2g 2 uo and is due to the inco- 
herent SP dynamics. The experimental data can be fitted reasonably well with the LP 
and SP spectra below 0.8 eV. Indeed we note that for high frequencies the effects due 

to the exchange-split bands and to the local Coulomb repulsion should be included.0 Thus 
the results can be considered meaningful for frequencies up to the absorption peak of the 
SP band. 



IV. EFFECTS OF AN EXTERNAL MAGNETIC FIELD 

In this section we discuss the effects of an external magnetic field on the optical properties 
of the system at x = 0.3. 

We first note that a unit of magnetic field ho can be fixed as 

= (30) 

If we choose for uq a reasonable value of 50meV, this unit is huge: in fact ho is of the 
order of 300 T. For realistic values of the magnetic fields the new energy scale g s fiBSho is 
small when compared with ooo- However in presence of an external magnetic field the phase 
diagram shows the tendency for LP — FM and FM — PI regions to grow. By introducing an 
external magnetic field (see Fig. 2), the magnetization does not vanish at T c , the transition 
temperature in absence of the field, showing the behavior of a first order transition.il'0 Thus, 
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in the vicinity of the transition temperature, the subtle balance between the two coexisting 
phases can be readily influenced by varying the magnetic field. 

The increase of the coherent motion due to the introduction of the external field is clearly 
shown in the calculated conductivity. We limit ourselves to the diagonal component. In Fig. 
3(a) we can notice the rise of the Drude term with increasing the magnetic field, while in Fig. 
3(b) the small-polaron absorption band is suppressed. This behavior is more pronounced 
at higher temperatures (see Fig. 4 for T = 0.96T C ). At T c the introduction of the magnetic 
field is able to fill up the low frequency optical gap present in the high-temperature phase 
(Fig. 5). Hence the magnetic field induces the transfer of spectral weight from high to low 
energies near the phase boundary.0 

The introduction of the magnetic field has dramatic consequences on the resistivity 
p = l/<r a ,a near T c , where a a ^ a results from the conductivity a a ,a(y}) in the limit uj — > (Fig. 
6). We obtain two different behaviors: metallic (dp/dT > 0) and insulating (dp/dT < 0) 
emphasizing in logarithmic scale the behavior of activated resistivity characteristic of a SP 
phase for T > T c .i0 The magnetic field shifts the peak of the resistivity at higher temper- 
atures. Since the PS regime is characterized by a rapid increase of the resistivity,^ at T c a 
large difference between resistivities at different magnetic fields occurs. The magnetoresis- 
tance ratio (see in Fig. 7(a) MR in percent) 

Mm = w^3_zm (31) 

p{ti) 

assumes large values at the transition temperature, so that the effect of the magnetic field 
is dramatically amplified near T c .0i Therefore, the CMR effect and, in general, the high 
sensitivity of the system to external parameters are due to the subtle balance between 
LP — FM and SP — PI phases characterized by different lattice, spin and charge properties. 

We have calculated the effect of the applied magnetic field on the DW factor. In Fig. 
7(b) we report A(DW), the relative variation of the DW factors obtained without and with 
the external field, in order to make a comparison with the magnetoresistance ratio. E1E1 The 
quantity A(DW) shows a well defined peak around the transition temperature, indicating 
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that the application of the magnetic field reduces the polaronic distortions. The two quan- 
tities, MR and A(DW) in Fig. 7, are both peaked around T c showing the existence of a 
strong correlation between lattice distortions, conductivity and magnetic properties.^ 



V. ISOTOPE EFFECT 

In this section we deal with the isotope effect in the CMR regime at x = 0.3. In 
particular we calculate the effects of isotope substitution on the infrared spectra. 

The isotopic substitution of the oxygen from 16 to l8 changes the values of Uq and g 
to ujI = ujq^M/M* and g* = g (M*/M) 1/4 respectively, where M* is the mass of 18 0. The 
coupling parameter A = g 2 uo/Qt is clearly unchanged. 

In Fig. 8 the phase diagram for l8 is presented and compared with that of 16 0. It is 
qualitatively altered with respect to the phase diagram of 16 showing a reduction of the 
FM regions .HHElHl! The heavier ion mass of 18 reduces the effective transfer integral t e ff 
(see eq. (|6])), so the region of coexisting phases and in particular the Curie temperature T c 
decrease. Indeed, in the limit where the intra-atomic exchange integral J is far larger than 
the bare transfer integral t, T c oc t e jj.i@ In the inset of Fig. 8, T*, the Curie temperature 
for 18 0, is reported along with T c , the Curie temperature for 16 0. We have evaluated the 
oxygen-isotope exponent a = —AlnT c / AlnM = — ( - Tc ~ Tc ^ ^ M ^_ M y making a comparison 
with experimental values.il We find the decreasing behavior of cto with increasing the hole 
doping x (Fig. 9). 

Near T c we have calculated the optical conductivity for the two isotopes (Fig. 10). While 
the 16 system is in the PS regime, 18 system is in the insulating phase. In the first case 
the optical response still shows the Drude term, for 18 only the small-polaron band at high 
frequency is present. Again this causes large changes in the conductivity of the system. 
Therefore the isotope effect represents another example of the sensitivity to the variation in 
external parameters. 
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VI. CONCLUSIONS 



We have discussed the effects of an external magnetic field and of the oxygen isotopic 
substitution on the infrared absorption spectra for x = 0.3 mainly near the phase boundary. 
Furthermore we have dealt with the lattice effects and the DW factor without and with the 
magnetic field. 

First we have focused our attention on LDPDF finding that, near and in the PI phase, 
only the sites occupied by an electron show strong lattice distortions and that the DW factor 
increases rapidly with temperature near T c in agreement with XAFS measurements.li!~il 

Then we have examined the consequences of an external magnetic field on the infrared 
absorption spectra at x = 0.3 stressing that, near the phase boundary, the application of 
the field induces a transfer of spectral weight from high to low energiesB and evaluating 
the CMR ratio. In our scheme, the found CMR is due to the subtle balance between 
coexisting phases characterized by different lattice, spin and electronic properties. The 
relation between the CMR ratio and the variation of the DW factor in presence of the 
magnetic field has pointed out the correlation that involves lattice distortions, conductivity 
and magnetic properties. 

Finally, within our approach, the oxygen isotope effect has been explained by the en- 
hanced electron-phonon coupling that induces a reduction of the FM phases. In the vicinity 
of the MI transition the isotope substitution induces large changes in the optical response 
and in the resistivity in agreement with experimental dataSilEllzl 

We believe that the strong sensibility of the infrared absorption with respect to the 
magnetic field and isotopic substitution could be used to test the validity of the PS scenario 
adopted in this paper. Indeed the experimental study of the infrared absorption is well 
established and it could be able to observe the effects discussed above. 

In our work, near the MI transition, the interplay of the electron-phonon interac- 
tion and the magnetic effects gives rise to a charge-density instability that takes into ac- 
count experimental evidences for coexisting localized insulating and delocalized metallic 
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component s.Elffl We have assumed that the scales of the inhomogeneities are much larger 
than the inter-particle distance and that the term arising from the mixing energy of the two 
coexisting phases and including the surface energy cost is not able to change qualitatively the 
behavior of the macroscopic quantities discussed. We have exploited a macroscopic PS sce- 
nario where the coexisting fractions of volume are determined by the variational procedure 
and the Maxwell construction. This approach to manganites provides results consistent with 
experimental measurements@@ and an explanation of the strong sensitivity of the system 
to external parameters. 

FIGURE CAPTIONS 

Fl (a) The distribution function of the atomic displacement (in units of ^ Mujq/2) at 
x = 0.3 as a function of the displacement coordinate X (in units of 2/^Mujq different 
from the length unit introduced above, the lattice constant a). 

(b) The Debye-Waller factor (in units of 2/Mujq) at x — 0.3 as a function of the 
temperature. 

F2 The magnetization (in units of saturation magnetization M$) as a function of the 
temperature for two different magnetic fields (in units of ho). 

F3 (a)-(b) The diagonal conductivity at T = 0.9 T c at different ranges of the frequency uj 
for several magnetic fields (in units of ho). The conductivities are expressed in units 
of e 2 c/mLdo, with c hole concentration and m = l/2t. 

F4 (a)-(b) The diagonal conductivity (in units of e 2 c/muo, with c hole concentration and 
m — l/2t) at T = 0.96 T c at different ranges of the frequency u for several magnetic 
fields (in units of ho). 

F5 The diagonal conductivity up to 18 ujq at T = 1.02 T c for different magnetic fields (in 
units of ho). The conductivity is expressed in units of e 2 c/muJo, with c hole concen- 
tration and m — l/2t. 
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F6 The inverse of the diagonal conductivity as a function of the temperature at different 
magnetic fields (in units of h ) (we have used u = 50 meV and the lattice constant 
a = 0.4 nm). 

F7 (a) The magnetoresistance ratio at a fixed magnetic field (in units of ho) as a function 
of the temperature. 

(b) The relative variation of the Debye- Waller factor at a fixed magnetic field (in units 
of ho) as a function of the temperature. 

F8 The phase diagram for the two different oxygen isotopes ( 16 dashed line, 18 O dot- 
ted line ). PI means Paramagnetic Insulator, FM Ferromagnetic Metal and AFI 
AntiFerromagnetic Insulator. The areas PI + FM and AFI + FM indicate regions 
where localized (PI or AFI) and delocalized (FM) phases coexist. In the inset the 
ferromagnetic transition temperatures as a function of the hole doping (T c dashed line 
for 16 and T* dotted line for 18 O). The temperatures are expressed in units of lu and 
the model parameters are t = 2uo, g = 2.8 (A = 0.65 for both isotopes) and e = O.Oli. 

F9 The oxygen-isotope exponent ao (circles) compared to the experimental exponent (di- 
amonds) deduced by Ref. 33 as a function of the hole doping. 

F10 The conductivity at T = 0.97 T c up to 19 lu for 16 and 18 0. The conductivities are 
expressed in units of e 2 c/mu , with c hole concentration and m — 1/2*. 
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